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Abstract 

We consider the problem of distortion minimal morphing of n- 
04 ' dimensional compact connected oriented smooth manifolds without 

, boundary embedded in K."+^. Distortion involves bending and stretch- 

ing. In this paper, minimal distortion (with respect to stretching) is 
, defined as the infinitesimal relative change in volume. The existence of 

1/^ ' minimal distortion diffeomorphisms between diffeomorphic manifolds 

, is proved. A definition of minimal distortion morphing between two 

isotopic manifolds is given, and the existence of minimal distortion 
, morphs between every pair of isotopic embedded manifolds is proved. 
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■ 1 Introduction 



A morph is a transformation between two shapes through a set of interme- 
diate shapes. A minimal morph is such a transformation that minimizes 
distortion. 

There are important applications of minimal morphing in manufactur- 
ing [71 [12] , computer graphics [THl [H] , movie making [9] , and mesh construc- 
tion [5l|6] . To address these applications would require a theory of minimal 
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morphing that includes bending and stretching together with algorithms to 
compute minimal morphs. In this paper, we do not address applications; 
rather, we formulate and solve the mathematical problem of minimal mor- 
phing with respect to stretching. A complete theory of minimal morphing 
is a subject for future research. 

In the research literature available at present, minimal morphing is con- 
sidered as a numerical problem where a cost functional is minimized over 
a finite number of intermediate shapes. We introduce a theory of distor- 
tion minimal morphing over a continuous family of states in the context of 
morphs between n-dimensional oriented compact connected smooth mani- 
folds without boundary embedded in M""^^ whose orientations are inherited 
from the usual orientation of M""*"^. The natural cost functional (for stretch- 
ing) measures the total relative change of volume with respect to a family 
of diffeomorphisms that defines the morph. This functional is invariant un- 
der compositions with volume preserving diffeomorphisms; hence, the cor- 
responding minimal morphs are not unique. On the other hand, we prove 
that the extremals of our functional are (strong) minima. Our main result 
is the existence of a distortion minimal morph (with respect to stretching) 
between every pair of isotopic submanifolds. 



2 Minimal distortion diffeomorphisms 

In this section we prove the existence of distortion minimal diffeomorphisms 
between diffeomorphic n-dimensional oriented manifolds M and (which 
are not necessarily embedded in R"+^) with respective volume forms lum 
and lon- 

Recall that the Jacobian of a diffeomorphism h : M ^ N is defined by 
the equation 

h*UJN = J{(^M,UJN){h) UJM, 

where h*LVN denotes the pullback of the volume form lon on A by the 
diffeomorphism h (see [l]). The Jacobian J(/i) := J{ujM,^N){h) depends 
on the diffeomorphism and the volume forms. 

The distortion (due to stretching) ^(m) at m G M, with respect to a 
diffeomorphism /i : M — > A, is defined by 

am) = H I^MAe) I I A, 1 ^ J(h){m) - 1, (1) 

where C M, for e > 0, is a nested family of (open) neighborhoods of the 
point m £ M such that Aa ^ whenever a > f3 > and H^^qA^ = m. 
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In other words, the distortion is the infinitesimal relative change of vol- 
ume with respect to h. It is easy to see that the definition of distortion does 
not depend on the family of nested sets Af,. 

We denote the set of all diffeomorphisms between manifolds M and 
N by Diff2(M, A/"). The total distortion functional $ : Diff^ (M,iV) R, 
with respect to the oriented manifolds (M, cjm) and {N,un), is defined by 

Hh) = J^(\j{h)im)\-iyiOM- (2) 

We will establish necessary and sufficient conditions for a diffeomorphism 
h : M N to he a minimum of the functional Also, we will show that 
a minimum always exists in Diff^(M, A^) provided that the manifolds are 
compact, connected, and without boundary. 

As a useful notation, we let X(M) denote the set of smooth vector 
fields on the manifold M. Also, we recall a basic fact from global non- 
linear analysis: Diff^(M, A^) is a Banach manifold and its tangent space 
at h e DiS^{M,N) can be identified with X{N). Indeed, an element of 
T/j Diff ^ (M, A^) is an equivalence class of curves [/ig], represented by a fam- 
ily of diffeomorphisms he with ho = h, where two curves passing through h 
are equivalent if they have the same derivative at h. For each n e N, this 
family defines a curve e <—>■ he{h~^ (n)) in N that passes through n at £ = 0; 
hence, it defines a vector Y G TnN by 

Y{n) := ^h,{h-\n)) 

The vector field Y G X{N) is thus associated with the equivalence class [h^]- 
In fact, the vector field Y does not depend on the choice of the representative 
of the equivalence class. On the other hand, for Y G X{N) with flow (j)t., we 
associate the curve ht = (pt ° h in Diff ^ (M, A'^) . The (tangent) equivalence 
class of this curve is an clement in ThBiS^{M, N). 

Proposition 2.1 (Euler-Lagrange Equation). Suppose that M and N are 
smooth connected compact orientable manifolds without boundary. A dif- 
feomorphism h : M ^ N is a critical point of the total distortion functional 
$ if and only if J{h) is constant. 

Proof. Let h^ : (—1, 1) T>iS'^{M, N) be a curve of diffeomorphisms from 
M to A^ such that ho = h. By definition, h G DiE^{M, N) is a critical point 

of the functional $(/i), if ^^(^t) = 0. Using the easily derived formula 
$(/i) = f J{hfujM - 2 Vol(iV) + Vol(M), 
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we note that /i is a critical point if and only if 



M 

Moreover, using the calculus of differential forms (see [Ij and note in particu- 
lar that Ly is used to denote the Lie derivative in the direction of the vector 
field y), we have that for = ipt o h, where ipt is the flow of y E X{N), 



d 
di 



_ d 
t=o dt 



{{ipt o h)*ujN) 



t=0 



t=Q 
t=0 



= h* LYU!j\f 

= h*(divYLON) 

= (divy) o hJ{h)ujM- 

By Stokes' theorem and the properties of the A-antiderivations d and iy-, 
we have that 



/ J{hf di-vY ohujM = I Jih) oh-^divYiON 
Jm Jn 

J{h) o LyojN = / J{h) o h~^diyu)N 



d{J{h)oh ^AiyuJ]\f)— / d{J{h)oh ^) AiyuJN 



N 



N 



N 



iy{d{J{h)oh ^)Aujn)+ I iy{d{J{h)oh ^))ujn 
d{J{h)oh-^){Y)uJN- 



N 



N 



Hence, h G Diff^(M, A^) is a critical point of the functional ^{h) if and only 
if 

d{J{h)oh-^){Y)ujN = Q 



JN 



for all y G X{N). It follows that if J{h) is constant, then /i is a critical 
point of <I>. 

To complete the proof it suffices to show that if 



/ df{Y)uJN = ^ 

Jn 



(3) 
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for all Y G X{N), then df = 0, where / := (J(/i) o /j-i). 

Suppose, on the contrary, that there exists a continuous vector field 
Y € X{N) such that df{Y){n) ^ for some point n ^ N. Without loss of 
generality, we assume the inequality df{Y)(n) > 0. The map df(Y) : N 
is continuous. Therefore, there exists an open neighborhood U C N oi the 
point n E so that df{Y){p) > for every p ^ U . Using the standard 
bump function argument (see PQ), we can construct a vector field Z E X{N) 
supported in U such that Jj^ df{Z)u)N = fjj df{Z)ui]\f > 0, in contradiction 
to equality ([3]). Hence, df = 0. □ 

Definition 2.2. A function h G Diff^(M, A^) is called a distortion minimal 
map if it is a critical point of the total distortion functional 

As an immediate corollary of proposition 12.11 we have the following the- 
orem. 

Theorem 2.3. A function h G DiE^{M,N) is a distortion minimal map if 
and only if J{h) is the constant function with value Vol(A^)/ Vol(M). 

We will use the elementary properties of distortion minimal maps stated 
in the following lemma. The proof is left to the reader. 

Lemma 2.4. Compositions and inverses of distortion minimal maps are 
distortion minimal maps. 

Also we will use (the strong form) of Moser's theorem on volume forms, 
which we state here for the convenience of the reader (see [8]). 

Theorem 2.5. Let Tt be a family of volume forms defined for t € [0, 1] on 
a compact manifold M. If 



for every n-cycle c on M , then there exists a one-parameter family of dif- 
feomorphisms cpt : M ^ M such that 



and 4>Q is the identity mapping. Moreover, the dependence of (ptinT-) on 
m (z M and t G [0, 1] is as smooth as in the family Tt. 

Theorem 2.6. If {Ad, com) o-nd [N,ojn) are diffeomorphic n-dimensional 
compact connected oriented manifolds without boundary, then (i) there is a 
distortion minimal map from M to N , (ii) every distortion minimal map 




(4) 



(t^ln = To 



(5) 
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from M to N minimizes the functional ^, and (Hi) the minimum value of 
^ is 

_ (Vol(M)-Vol(iV))^ 
Vol(M) ■ 

Proof. To prove (i), choose a diffeomorphism h £ Diff^(M, N) and note that 
the differential form h*uJN is a volume on M. Define a new volume on M 
as follows: 

Vol(M^ 
Jm ^ 



Since 



/ = I 

Jm Jm 



and M is compact, by an application of Moser's theorem 12.51 there exists a 
diffeomorphism f : M —> M such that ujm = f*^M- Hence, 

JMh*UJN , 
-y^^M = {hof) 

and, since M is connected, we conclude that J{h o /) = Vol(A^)/ Vol(M) is 
constant and DJ{h o /) = 0. Thus, k = ho f \s a, distortion minimal map. 

To prove parts (ii) and (iii), note that if k is an arbitrary distortion 
minimal map from M to A^, then 

2^r.,..-s (Vol(M) - Vol(7V))^ 



^k) = {\J{k) \ - l)^Vol(M) 



Vol(M) 



We claim that this value of $ is its minimum. 

Let g G Diff^(M, A^). By the Cauchy- Schwartz inequality, 

m= I {\J{g)\-ifu:M 

JM 

1 X \2 



voi(M) V y^,^ 

^ (Vol(M) - Vol(iV))^ 



Vol(M) 



as required. □ 
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Example 2.7. Let Sr and Sji be two-dimensional round spheres of radii r 
and R (respectively) centered at the origin in M."^. Define h : Sr ^ Sr hy 
h{p) = R/rp for p = (x,y,z) G 5*^. We will show that h is a distortion 
minimal map. 

Let ujr (respectively, wr) be the standard volume forms on Sr (respec- 
tively, Sr) generated by the usual Riemannian metric gp{X,Y) = {X,Y) 
for X,Y e 

Using the parametrization of Sr and Sr by spherical coordinates, it is 
easy to compute that the Jacobian 



for all m G Sr', hence, by theorem 12.61 ^ is a distortion minimal map. 

Remark 2.8 (Harmonic maps). For h € Diff^(M, A^), the distortion func- 
tional ([2]) has value 



Thus, it suffices to consider the minimization problem for the reduced func- 
tional ^ given by 



We note that if M and are one-dimensional, then ^ is the same as 



An extremal of this functional is called a harmonic map (see OEISI). Thus, 
for the one-dimensional case, distortion minimal maps and harmonic maps 
coincide. On the other hand, there seems to be no obvious relationship in 
the general case. 

3 Morphs of embedded manifolds 

We will discuss a minimization problem for morphs of compact connected 
boundaryless oriented n-dimensional smooth manifolds embedded in M"+^. 



J{uJr,ujR){h){m) = R^/r^ = \o\{Sr) /Yo\{Sr) 



^{h) = 1 \J{h)\^UM - 2Vol(iV) + Vol(M) 



JM 
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3.1 Pairwise minimal morphs 

Definition 3.1. Let M and N be compact connected oriented n-dimensional 
smooth manifolds without boundary embedded in M"^^. A function 
if : [0, 1] X M ^ W^^^ is a morph from M to iV if the following conditions 
hold: 

(i) p H{t,p) IS a diffeomorphism onto its image for each t G / = [0, 1]; 

(ii) the image M* = H(t, M) is an n-dimensional manifold possessing all 
the properties of M and iV mentioned above; 

(iii) p I— > H{0,p) is a diffeomorphism of M; 

(iv) the image of the map p H{l,p) is A^. 

We denote the set of all morphs between manifolds M and N by A4{M, N). 

For simplicity, we will consider only morphs H such that p H{0,p) 
is the identity map. We assume that each manifold M* = H{t, M) (with 
= M and = N) is equipped with the volume form ujt = irn^, where 



is the standard volume form on and rjt : M* M""*"^ is the outer unit 

normal vector field on M* with respect to the usual metric on R""*"^. Also, 
as a convenient notation, we use /i* = H{t, •) : M — > M*. 

Definition 3.2. A morph H is distortion pairwise minimal (or, for brevity, 
pairwise minimal) if /i*'* = /i* o {h^)~^ : M* is a distortion minimal 

map for every s and t. We denote the set of all distortion pairwise minimal 
morphs between manifolds M and N by 'PM{AI, N). 

By proposition 12.11 and theorem 12.61 ^ morph H is pairwise minimal if 
and only if each Jacobian J{uJs,uJt){h^'^) is constant. 

Proposition 3.3. Let M = and N = AI^ be n-dimensional manifolds 
as in definition \3.1\ equipped with the (respective) volume forms ujq and loi. 
A morph H between M and N is distortion pairwise minimal if and only if 



n = dxi A dx2 A ... A dxn-^-i 



J{uQ,ujt){h*){m) 
Vol(M*) 



Vol(M) 



1 



(7) 



for all t G [0, 1] and m £ M. 
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Proof. Using lemma 12.41 and theorem 12.6^ it suffices to prove that each map 
: M ^ M* is minimal if and only if the map d?]) is constant. An applica- 
tion of theorem 12.31 finishes the proof. □ 

Proposition 3.4. Let M and N be n- dimensional manifolds as in propo- 
sition \3.3l If there is a morph G from M to N , then there is a distortion 
pairwise minimal morph between M and N . 

Proof. Fix a morph G from M to with the corresponding family of dif- 
feomorphisms 5* := G{t,-), let M* := G{t,M), and consider the family of 
volume forms 

Vol(M) , , 
'"* = Vol(M^^^)^* 

defined for t € [0, 1]. It is easy to see that 



= / Wo; 
M Jm 

hence, by Moser's theorem 12. 5| there is a family of difFeomorphisms a* on 
M such that a* depends continuously on t and ujm = (a*)*Wi. It follows 
that 

. , . , Vol(M*) 

therefore, 

Vol(M*) 



J{uJM,i^t){g o a )(m) 



Vol(M) 

for all m € M. The morph H corresponding to the family /i* := (7* o a* is 
the desired distortion pairwise minimal morph. □ 



3.2 Minimal morphs 

We will define distortion minimal morphs between embedded connected ori- 
ented n-dimensional smooth manifolds without boundary. 

For a morph H from M to N, let Es^t denote the total distortion of 

Es,t= [ (\j{h''')\-lYuJs 
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If is a morph, then Es^t is twice continuously differentiable with respect 
to s. By Taylor's theorem, Eg^t has the representation 

Es,t = E,, + ^^{E^Aa^ - + ^^(^«-*)L=*(^ - + 0{is- tf). 
Note that E^^t and both vanish, and 

Definition 3.5. The infinitesimal distortion of a morph from M to 
TV at t G [0, 1] is 

\ 2 



E., r mh^) 



s^t {s - tf A 



The total distortion functional ^ defined on such morphs is given by 

2 



Definition 3.6. A morph is called a distortion minimal extremal if it is 
an extremal of the functional ^> with respect to morphs. A morph is 
called a distortion minimal morph if it minimizes the functional ^. 

Note that £^ (t) depends continuously on t provided that H is a 
morph. 

Lemma 3.7. For every morph H G JV[[M,N) there exists a pairwise min- 
imal morph G G VM{M,N) such that < ^{H). In particular, if 
H G A4{M,N) is a distortion minimal morph, then there exists a pair- 
wise minimal morph G G 'PA4{M, N) such that ^{H) = ^(G). 

Proof. Let H G Ai{M,N) be a morph with the intermediate states M* = 
H{t,M). By proposition 13.41 there exists a pairwise minimal morph G G 
VA4{M, N) with the same intermediate states. The deformation energy of 
transition maps satisfies the inequality Es^t{H) > Es^t{G) for all s,i G [0, 1] 
because G is pairwise minimal. Therefore, £^ {t) > for all t G [0,1], 

and, consequently, 

^{H) > (9) 
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as required. 

If H is distortion minimal, the inequality ^{H) < ^{G) holds. Com- 
paring the latter inequality with inequality Q, we conclude that ^{H) = 

Corollary 3.8. (i) The following inequality holds: 

inf $(G) < inf ^(H). (10) 

GeVM{M,N) H€M{M,N) 

(a) If there exists a minimum F of the total distortion functional ^ over 
the class 'PA4{M, N), then F minimizes the functional <1> over the class 
M{M,N) as well: 

= min ^(G) = min <^(H). (11) 

GeVM{M,N) H€MiM,N) 

Lemma 3.9. The total distortion of a pairwise minimal morph H from 
M to N is 

Proof. The proof is an immediate consequence of formula ([8]) and proposi- 
tion [331 □ 



Lemma 3.10. Consider an auxiliary functional 

*(</.) = ^dt (13) 

defined on the admissible set 

Q = G Ci([0, 1];M+) : m = \ol{M) , <f>{l) = Vol(iV)}. 

The functional (|13p attains its minimum at(j) = Vol(M) i/Vol(M) = Vol(A^) 
and at 4){t) = (y^VoI(M)- Y^Vol(iV) )t-^^\o\{M) whenever the volumes 
of M and N are distinct. The minimal value of the functional ^ is 

min^'(/9) = 4(VVol(iV) - VVol(M))^ (14) 

Proof. The proof is a simple application of the Euler-Lagrange equation and 
the Cauchy-Schwarz inequality. □ 
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Using corollary 13.81 and lemma [3.101 we will minimize the total distortion 
energy functional <I> over the set A4{A4, N) of all morphs. 

Theorem 3.11. Let M and N be two n- dimensional manifolds satisfying 
the assumptions of definition \3. 1[ If M and N are connected by a morph, 
then there exists a distortion minimal morph. The minimal value of $ is 

min ^(H) = A(J\o\(N) - J\o\iM)f . (15) 
H&M{M,N) V V V / V V // 

Proof. Let G be a morph between M and A^. Without loss of generality, we 
assume that G is pairwise minimal (see proposition 13. 4p . Set 

H{t,m) = X{t)G{t,m), 

where A : [0, 1] ^ M is to be determined. 

Note that if M* = H{t,M) and = G{t,M), then 

Vol(M*)= / (/i*r^M= [A(t)]" / (5*)*a;M= [A(t)]"Vol(VF*)- 
Jm Jm 

Let (p{t) be the minimizer of the auxiliary functional ^ from lemma 13.101 
and define 

X(t) - 

^ ^ " Lvol(H^*) 

The corresponding volume Vol(M*) = (p{t); therefore, by corollarv 13.81 and 
lemma [3.10( the morph H minimizes the total distortion functional $ over 
the class M{M,N) and ^{H) = 4(v^Vol(iV) - ^Vol(M))^ □ 

The next result provides a basic class of distortion minimal morphs. 

Proposition 3.12. Suppose that M is an n-dimensional manifold embedded 



m 



that satisfies the assumptions of definition \3.1[ If a is a positive 



real number and 

N := {am : m £ M}, 

then 

(i) N is a manifold satisfying all of the assumptions of definition \3.1 
(a) the morph given by the family of maps h^{m) = X{t)m, where 



X{t) = Yol{M)-n (/Vol(M) - ^Vol(7V) )t - VVol(M) 
is distortion minimal. 
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Proof. Define /i*(m) = X{t)m. It is easy to check that /i* defines a morph 
from M to N. Also, we have that J(/i*) := J{oJM,^t){h*) = [-^(i)]"- Since 
J(/i*) is constant on M, the family /i* defines a pairwise minimal morph H. 

We will determine A(t) so that the morph H becomes a minimizer of ^ 
over the class ^A{M,N). Indeed, by lemma lB.lUl it suffices to choose A so 
that 

Vol(M*) = [A(t)]"Vol(M) = [(\/Vol(M) - ^yYol{N) )t - VVol(M) 
which yields 



A(t) = Vol(M)-- ( VVol(M) - ^/Yo\{N) )t - ^/Yo\{M) 

The corresponding morph H{t, m) = \{t)m satisfies the equality 

^(H) = min ^>(G). 

G(^M{M,N) 



□ 
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